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We present a planar qubit design based on a superconducting circuit that we call concentric transmon. While 
employing a straightforward fabrication process using A1 evaporation and lift-off lithography, we observe 
qubit lifetimes and coherence times in the order of 10 /iS. We systematically characterize loss channels such as 
incoherent dielectric loss, Purcell decay and radiative losses. The implementation of a gradiometric SQUID 
loop allows for a fast tuning of the qubit transition frequency and therefore for full tomographic control of 
the quantum circuit. Due to the large loop size, the presented qubit architecture features a strongly increased 
magnetic dipole moment as compared to conventional transmon designs. This renders the concentric transmon 
a promising candidate to establish a site-selective passive direct Z coupling between neighboring qubits, being 
a pending quest in the field of quantum simulation. 


Quantum bits based on superconducting circuits are 
leading candidates for constituting the basic building 
block of a prospective quantum computer. A common 
element of all superconducting qubits is the Josephson 
junction. The nonlinearity of Josephson junctions gen¬ 
erates an anharmonic energy spectrum in which the two 
lowest energy states can be used as the computational 
basis^’^. Over the last decade there has been a two order 
of magnitude increase in coherence times of supercon¬ 
ducting qubits. This tremendous improvement allowed 
for demonstration of several major milestones in the pur¬ 
suit of scalable quantum computing, such as the con¬ 
trol and entanglement of multiple qubits^’^. Further in¬ 
creases in coherence times will eventually allow for build¬ 
ing a fault tolerant quantum computer with a reasonable 
overhead in terms of error correction, as well as imple¬ 
menting novel quantum simulation schemes by accessing 
wider experimental parameter ranges^. While supercon¬ 
ducting qubits embedded in a 3D cavity® have shown 
coherence times in excess of 100/xs^, this approach may 
impose some constraint on the scalability of quantum cir¬ 
cuits. Since the Josephson junction itself does not limit 
qubit coherence®, comparably long lifetimes can also be 
achieved in a planar geometry by careful circuit engineer¬ 
ing. 

In this paper, we present the design and characteri¬ 
zation of a superconducting quantum circuit comprising 
a concentric transmon qubit®, schematically depicted in 
Fig. 1(a). The two capacitor pads forming the transmon’s 
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large shunt capacitance are implemented by a central 
disk island and a concentrically surrounding ring. The 
two islands are interconnected by two Josephson junc¬ 
tions forming a gradiometric SQUID. A SOLI impedance 
matched on-chip flux bias line located next to the qubit 
allows for fast flux tuning of the qubit frequency due to 
the imposed asymmetry. This guarantees high experi¬ 
mental flexibility and enables full tomographic control. 
The gradiometric flux loop design reduces the sensitivity 
to external uniform magnetic fields and thus to external 
flux noise. For readout and control purposes we embed 
the qubit in a microstrip resonator circuit, forming a fa¬ 
miliar circuit quantum electrodynamics (cQED) system®. 
The sample is fabricated in hybrid coplanar-microstrip 
geometry, featuring a ground plane on the backside of 
the substrate, thus exploiting the increase in mode vol¬ 
ume in a microstrip architecture. The circular shape of 
the transmon features a strongly reduced electric dipole 
moment due to symmetry and hence a reduction in radi¬ 
ation loss can be expected^® compared to a regular pad 
geometry of equal size and electrode spacing. The sample 
is prepared in a straightforward fabrication process using 
pure aluminum for the metalization and employing the 
conventional shadow angle evaporation technique. The 
best measured relaxation and decoherence times are in 
the order of 10 fas. 

The schematic circuit diagram of the concentric trans¬ 
mon is depicted in Fig. 1(c). The two Josephson junc¬ 
tions are connecting the central island to the outer ring 
of the transmon. The islands act as coplanar electrodes^® 
giving rise to the total qubit capacitance C = 81 fF, in¬ 
cluding the contribution by the ground plane and cou- 
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FIG. 1. (a) Geometry of the concentric transmon qubit. Two 
Josephson junctions (red crosses), located opposite to each 
other, connect the central island to an outer ring, (b) Optical 
micrograph of the fabricated sample. The readout resonator 
(red) capacitively couples to the concentric transmon from 
above. The on-chip flux bias line (blue) is visible to the right. 
It is designed in coplanar geometry, having the microwave 
ground reference in the device plane (bright color). The flux 
bias line is grounded at one end on the chip, (c) Schematic 
circuit diagram of the concentric transmon, revealing the gra- 
diometric SQUID architecture. The central part, marked with 
a black dot, corresponds to the center island of the transmon. 
Since the mutual inductances to the flux bias line are not 
equal. Mi 7 ^ M 2 , the effective critical current of the SQUID 
can be tuned. 


pling capacitances. Considering loop Li as the pri¬ 
mary transmon loop, the gradiometric dc-SQUID archi¬ 
tecture can be recognized. The kinetic inductance of the 
superconducting aluminum wire can be neglected due 
to its thickness and the width of the Josephson junc¬ 
tion leads. The effective critical current Ic,eff of the 
SQUID is tuned by applying an inhomogeneous mag¬ 
netic field supplied by the on-chip flux bias line. Due 
to the gradiometric geometry, the effective critical cur¬ 
rent Ic,eff = 2/c| cos( 7 rA$/ 2 <i)o)| in the primary loop Li 
is 2$o-periodic in the flux asymmetry A$ = |$i—$ 2 ! be¬ 
tween the loops. Here, $0 = ^/2e denotes the magnetic 
flux quantum. We analytically calculate the net mutual 
inductance to the flux bias line to be 2.3 pH. This yields 
a flux bias current of 0.9 mA required to induce a <i>o in 
the primary transmon loop. 

The fabrication process consists of two subsequent de¬ 
position and lithography steps. The feedline, the mi¬ 
crostrip resonator and the flux bias line are structured 
by optical lithography in a lift-off process of a 50 nm 
thick aluminum film. After that, the concentric trans¬ 
mon including the Josephson junctions are patterned us¬ 
ing electron beam lithography. The Josephson junctions 
are formed by shadow angle evaporation with electrode 
film thicknesses of 30 nm and 50 nm, respectively. A 
50 nm thick aluminum film is applied on the backside 
of the intrinsic silicon substrate. Both metalizations on 
the device side of the chip are deposited by aluminum 
evaporation at a background chamber pressure of about 
3T0“® mbar. A micrograph of the fabricated sample is 


TABLE I. Calculated loss contributions fo the concentric 
transmon qubit. The main contribution arises from 

inductive coupling to the flux bias line, leading to a total 
Purcell limitation of rj"p = 16 ps. The estimated reciprocal 
sum F))^ is in good agreement with the measured Ti. 


Purcell 

defects 

radiation 

reciprocal sum 

p-i r~^ 

l,sm ^ l,ind ^ l,cap 

p-i 

^ 1,TLF 

p-i 
^ l,rad 

p-i 

^ s 

47 fis 32 ~ 87 /^s 




r“), = 16 /xs 

~ 26 /is 

> 100 tis 

8.9 fis 


depicted in Fig. 1(b). 

We operate the device at a fundamental qubit tran¬ 
sition frequency lOql'l'K = 6.85 GHz, away from any 
flux sweet spot and far detuned from the readout res¬ 
onator to reduce Purcell dissipation. The Josephson 
energy Ej jh = 29 GHz dominates the charging energy 
Ec/h = 0.24 GHz, yielding Ej/Eq = 120. This assures 
the circuit to be inherently insensitive to charge noise 
and features an intrinsic self-biasing^^. The total criti¬ 
cal current of both Josephson junctions is 58 nA. Due 
to its large loop size, the concentric transmon exhibits 
a notable geometric inductance Lg = 0.6 nH, accounting 
for 10% of the total qubit inductance. As one impor¬ 
tant consequence, the qubit anharmonicity is not con¬ 
stant with respect to the bias point. For the operation 
point we extract a relative anharmonicity of 3.4%. Rele¬ 
vant qubit parameters are extracted by fitting an effective 
system Hamiltonian to spectroscopic data, see Supple¬ 
mental Material^^. By sweeping the flux bias current we 
observe a period of 1.7 mA in the qubit transition fre¬ 
quency, confirming the expected 2 $o-periodicity in good 
approximation. 

Since coherence is relaxation limited, T 2 < 2Ti, tak¬ 
ing effort in increasing Ti typically results in an increase 
in T 2 . To accomplish this, we engineered the concentric 
transmon to have a reduced sensitivity to its major loss 
channels, namely spontaneous Purcell emission, radiative 
dipole decay and loss due to surface defect states. Losses 
due to quasiparticle tunneling processes typically impose 
a Ti limitation at around ~ Ims^^, having no consider¬ 
able effect on the lifetime of our circuit. Relevant contri¬ 
butions leading to qubit decay are summarized in Tab. I. 

The coupling limited quality factor Ql = of 

the dispersive readout resonator at Wr/27r = 8.79 GHz 
is 2.1-10^, close to the design value. The qubit lifetime is 
potentially Purcell limited by spontaneous emission into 
modes that are nearby in frequency. Major contributions 
are the dispersive single-mode decay into the capacitively 
coupled readout resonator as well as emission into the 
flux bias line^^d4 inductive coupling. The coupling 

constant g/27r = 55 MHz between qubit and resonator 
is extracted from the dispersive shift of the resonator^®. 
We calculate a single-mode Purcell limitation induced by 
the readout resonator of Pi^sm = = 1/47ps at 

a detuning A = 1.94 GHz from the readout resonator. 
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The multi-mode Purcell limitation due to inductive cou¬ 
pling of the qubit to the flux bias line of impedance 
Zq = son is Ti^ind = ujqM'^KLZo) = l/32^s, with 
M = 2.3 pH being the gradiometric mutual inductance 
and L = 6.3 nH the total qubit inductance. This rather 
stringent loss channel can be reduced further by decreas¬ 
ing the coupling to the flux bias line. By an analogous 
approach to account for Purcell decay due to capacitive 
coupling to the flux bias line {Cc ~ 0.1 fF), we calculate 
ri,cap = ojqC'^Zo/C ^1/87/xs. The presented values are 
approximations and give a rough estimate of the expected 
Purcell loss. A more stringent analysis would require a 
full 3D electromagnetic simulation of the circuit and a 
blackbox quantization^®. 

The coplanar concentric transmon is embedded in a 
microstrip geometry, where the ground reference on the 
device side of the substrate is substituted by a back¬ 
side metalization. The largest fraction of the electric 
held energy is stored in the substrate, and the held 
strength at incoherent and weakly coupled defects re¬ 
siding in surface and interface oxides of the sample is 
reduced due to an increased mode volume. Highest helds 
in the geometry appear in the gap between the center 
island and the ring of the concentric transmon within 
the substrate. From the vacuum energy of the trans¬ 
mon we extract a weighted mean electrical held strength 
|if| = 2.3 V/m in the surface and interface oxide of an 
estimated ehective volume oi V = 50T0“^®m®. We as¬ 
sume a maximum defect dipole moment |(io| = 1.6 eA, 
reported in literature^^’^® as the highest dipole moment 
observed in Josephson junction barriers and therefore 
yielding a worst case estimation. We employ a normal¬ 
ized dipole moment distribution^^ P{p) = A\/l— 
with relative dipole moment p = |d|/|(io|. Taking into 
account a normalized defect probability distribution^® 
P{uj,d) = cos°‘ 6/ s\x\9 one can estimate the mean 
relaxation rate Ti tlf due to a single incoherent two 
level huctuator (TLF) with averaged parameters to be 

fdo 2\p\2 riiiTLF /■’r /2 

j dpP{p)^^ j da; y d0PiLo,e)Cicop. 

( 1 ) 

to denotes the TLF frequency that we integrate to a maxi¬ 
mum frequency = 15 GHz, and 9 sets the dipole 

matrix element sin0. The spectral density C{ujq) = 
sin^ 072 /( 7 |-l-(a;—ojg)®) essentially is the Fourier trans¬ 
form of the coupling correlation function^®, with a TLF 
dephasing rate 72/271 = 10 MHz. The averaged rate in¬ 
duced by a single TLF, given in Eq. (1), is multiplied by 
the number N = poVHujtlf of defect TLF interacting 
with the qubit. With a distribution parameter a = 0.3®® 
and a constructed defect density po = 4-10®/pm®/GHz 
we compute a Ti limitation due to the bath of incoher¬ 
ent TLF to be Fi .tlf 1/26/xs. The choice of po is 
consistent with literature^^’®^ and is justified by a very 
good agreement with a loss participation ratio analysis 
carried out via a finite element simulation, see Supple¬ 
mental Material^®. The calculated decay rate shows a 



FIG. 2. (a) Measured relaxation time Ti = 9.1/xs and de¬ 
phasing time T 2 = 10 /xs with applied measurement pulse 
sequences shown in the inset. The operation frequency is 
a)q/27r = 6.85 GHz, corresponding to a detuning of 1.9 GHz 
below the readout resonator frequency. A Hahn echo pulse 
is applied in the Ramsey sequence to compensate for low- 
frequency fluctuations in the qubit transition frequency, (b), 
(c) Ti decay and bare dephasing without echo pulse T 2 mea¬ 
sured at the flux sweet spot at 7.72 GHz. While Ti is reduced 
due to Purcell dissipation, TJ is comparable to the bare Tp 
measured far detuned from the sweet spot. The slight devi¬ 
ation from an exponential decay is attributed to excitation 
leakage. 


very weak dependence on the employed parameters 72 , 
a and integration cutoffs. Fi tlf imposes a limitation 
for the best measured Ti due to the quasi-static bath of 
incoherent defects. In general, defects in the Josephson 
junction of the device do not significantly contribute to 
qubit decay. Due to its small size, the defect density in 
the Josephson barrier is discretized and therefore highly 
reduced®®’®®. 

As pointed out in Ref. [10], radiative loss becomes ap¬ 
parent for qubits with a large electric dipole moment. 
Radiative decay is reduced as the dipole of the mirror 
image, induced by the ground plane of the microstrip 
geometry, radiates in antiphase, leading to destructive 
interference. Our circular geometry brings about an ad¬ 
ditional decrease in radiated power by strongly reduc¬ 
ing the electric dipole moment of the qubit. We ana¬ 
lyze this by simulating the dissipated power in a con¬ 
ductive plane placed 1.5 mm above our geometry in the 
medium far field and compare the result to a conven¬ 
tional pad architecture. The internal quality factor of 
the qubit eigenmode indicates a radiative contribution of 
^irad ^ 100/xs. The radiative dissipation of a compa¬ 
rable qubit in pad geometry exceeds this value by more 
than an order of magnitude. 

The dissipative dynamics of the investigated concen¬ 
tric transmon is depicted in Fig. 2. We excite the qubit 
by applying a tt pulse and measure its state (dz) with a 
strong readout pulse after a varying time delay At. The 
obtained energy relaxation time Ti = 9.1 /xs in Fig. 2(a), 
lower part, is in good agreement with the estimated Ti 
time, see Tab. 1. To verify the presented loss participa¬ 
tion ratio analysis, Ti decay close to the flux sweet spot 
is shown in Fig. 2(b). The Purcell contribution is cal- 
























4 


culated to be 4.7 ^s, reducing the overall dissipation es¬ 
timate to 3.7 /iS. The anharmonicity at the flux sweet 
spot is reduced to 1 %, favoring leakage into higher levels 
and therefore degrading the reliability of the exponential 
fit. The reduction in Ti at the sweet spot is not fully 
understood. A possible explanation is that the presented 
first order model for estimating Purcell decay is not very 
accurate for this high geometric inductance transmon. 
In addition, losses related to the geometric inductance 
such as quasi particles or magnetic vortices may be more 
prominent at the flux sweet spot. 

In continuous lifetime measurements, we observe tem¬ 
poral fluctuations in Ti down to about 2 ^s. We conjec¬ 
ture that this is attributed to discrete TLF dynamics of 
individual TLF, located in the small oxide volume at the 
leads to the Josephson junctions, where the electric field 
strength is enhanced. 

The spin echo dephasing time is T 2 = 10 /xs, see 
Fig. 2(a), upper part. Due to the Hahn echo pulse in be¬ 
tween the projecting 7r/2 pulses, our device is insensitive 
to low-frequency noise roughly below 1/At = 25 kHz. A 
Ramsey T 2 ^ 2 /xs is measured without echo pulse. It is 
interesting to note that a rather high T 2 is maintained in 
spite of the large loop size and while operating the trans¬ 
mon detuned by 1 GHz from its flux sweet spot, where 
flux noise contributes to pure dephasing. Figure 2(c) 
shows the dephasing time T 2 without echo pulse mea¬ 
sured at the flux sweet spot. Dephasing, presumably 
induced in part by local magnetic fluctuators, is reduced 
due to the vanishing slope of the energy dispersion at 
the flux sweet spot. In a non-tunable TiN version of 
the concentric transmon with a single Josephson junc¬ 
tion and an opening in the outer ring, we find maximal 
Ti ~ 50 /xs and dephasing times up to 60 /xs without 
Hahn echo. Similar coherence times might be achieved 
for the tunable concentric design by implementing it in 
a TiN material system. The measured coherence times 
compare with other planar transmon geometries such as 
a non-tunable A1 based transmon with decreased finger 
gap size^"^ (Ti = 9.7 /xs, T| = 10 /xs) and the cross shaped 
transmon^^ (Ti = 40/xs, T 2 = 20/ts). 

Figure 3 demonstrates fast frequency control of the 
concentric transmon. This is commonly referred to as 
Z control since Pauli’s spin operator appears in the 
Hamiltonian representation of the qubit. We record the 
equatorial precession (see Fig. 3(b)) due to a detuning 
pulse of amplitude rj in between two projective 7r/2 
pulses. The pulse sequence is given in Fig. 3(a). The 
pulse amplitude 77 translates into magnetic flux applied 
to the flux bias line, see the experimental setup in the 
Supplemental Material^^. In the laboratory frame, this 
corresponds to a shift in qubit frequency A/ oc 77 , which 
is confirmed in a quasi-spectroscopic measurement by ex¬ 
citing the detuned qubit with a tt pulse, see Fig. 3(c). In 
the frame rotating with the qubit frequency the pulse 
induces a change in Larmor frequency, leading to an ef¬ 
fective precession in this rotating frame with an angular 
frequency proportional to the detuning pulse amplitude 



/ (GHz) 


FIG. 3. Demonstration of fast Z tunability of the concentric 
transmon. (a) Pulse sequence applied for calibration. Be¬ 
tween two projecting 7r/2-pulses we apply a Z-rotation of 
amplitude 77 and length At. (b) The expectation value of the 
qubit state oscillates between its fundamental states | 0 ), | 1 ) 
with a relative Larmor frequency ujl /27r = 65 MHz for a pulse 
amplitude of 77 = 32 /xA. (c) Check measurement of the bare 
and detuned qubit transition frequency tu, by exciting with a 
TT-pulse (see inset). The initial frequency is shifted dependent 
on the amplitude 77 of the applied current pulse, see different 
Lorentzian fits. For 77 = 32/xA, as applied in (a), we obtain a 
frequency shift in accordance with the one extracted in (b), 
see red Lorentzian fit. 


77 . By increasing 77 , we demonstrated a fast frequency 
detuning of up to 200 MHz (not shown here). A further 
increase of the tunability range, requiring pulse shaping, 
renders our device a valuable tool for a variety of quan¬ 
tum experiments. 

We propose the concentric transmon as a suitable can¬ 
didate to establish a direct inductive Z coupling be¬ 
tween neighboring qubits. Since the area of the mag¬ 
netic flux loop is large, the magnetic dipole moment of 
our qubit is strongly enhanced as compared to conven¬ 
tional transmon designs where it is typically negligible. 
For two concentric transmon qubits separated by 150 /xm, 
we analytically calculate a maximum mutual inductance 
Mij = 31 pH by applying the double integral Neumann 
formula^®. Comparing the ratio Mij/L ^ 0.5% to the 
ratio /? = Cg/C ~ 10%, typical for cQED circuits and 
corresponding to a 100 MHz coupling, we estimate an 
inductive Z coupling up to (7™‘^/27r ^ 5 MHz. Here, 
L = 6.3 nH is the total qubit inductance and /3 denotes 
the ratio of coupling capacitance Cg to total qubit ca¬ 
pacitance C. At the given qubit distance, the capacitive 
coupling is simulated to be in the order of 40 MHz and 
therefore less than a factor of ten higher compared to the 
inductive coupling. To further increase the Z coupling, 
the qubit geometry can be adapted to enhance the mu¬ 
tual inductance. Especially for neighboring qubits being 
detuned in frequency, the Z coupling may dominate the 
effective transversal capacitive coupling. Since the mu- 
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tual inductance vanishes completely for adjacent qubits 
arranged at a relative rotation angle of 7r/2 and is re¬ 
duced by more than an order of magnitude when their 
Josephson junctions are aligned, the concentric transmon 
geometry allows for a site-selective Z coupling. We con¬ 
sider this scheme as highly promising for the field of quan¬ 
tum simulation, for instance in the context of implement¬ 
ing Ising spin models. Patterning an array of concen¬ 
tric transmon qubits featuring X and Z coupling along 
two orthogonal directions while exploiting the strongly 
reduced off-site inductive interaction may be a route to 
implement a quantum neural network, which is a power¬ 
ful tool in quantum computation^®. 

In conclusion, we introduced a planar tunable qubit 
design based on a superconducting circuit that we call 
concentric transmon. The observed qubit lifetimes and 
coherence times are in the order of 10 /iS and thereby 
competitive with conventional transmon geometries. The 
qubit lifetime is Purcell limited by its readout resonator 
and the flux bias line. Radiative loss, which is an intrinsic 
loss channel to the geometry, is demonstrated to be limit¬ 
ing only above ~ 100 /rs, evincing the potential of the re¬ 
ported architecture. A major advantage of our approach 
is the straightforward fabrication process. We demon¬ 
strated full tomographic control of our quantum circuit 
and discuss the high potential of the presented qubit de¬ 
sign for the implementation of a direct site-selective Z 
coupling between neighboring qubits, being a pending 
quest in quantum simulation. 
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I. EXPERIMENTAL SETUP 


qubit 
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FIG. SI. Schematic diagram of the experimental setup. 


for additional magnetic shielding. The base stage of the 
cryostat is surrounded by a radiation shield having its 
inner surface coated with blackbody-absorbing material 
to shield the sample from stray infrared radiation gener¬ 
ating quasiparticles®^. 

The schematic diagram of the measurement and mi¬ 
crowave setup is depicted in Fig. SI. To perform mea¬ 
surements on the concentric transmon qubit we use a 
pulsed readout tone. The resonator shift is extracted 
from the microwave reflection signal at a single-ended 
50 n matched transmission line that capacitively couples 
to the readout resonator. We down-convert the read¬ 
out pulse by heterodyne sideband mixing to eliminate 
dc noise in the Fourier spectrum yielding the dispersive 
resonator shift. Qubit manipulation is done by hetero¬ 
dyne single sideband mixing. This guarantees that the 
free evolution of the qubit is not disturbed by microwave 
leakage through the IQ mixer during pulse off time. In 
that way other sidebands and the local oscillator (LO) 
leakage are reduced by at least 40 dB in power. For XY 
control of the qubit, microwave pulses are applied to the 
same transmission line used for the readout. 

The qubit transition frequency is adjusted by a DC 
current applied to the on-chip flux coil. High frequency 
noise is hltered at the 4 K stage with RCR type 7r-hlters 
at about 25 kHz and on the base plate via an RC-element 
enclosed in copper powder®^. Fast flux pulses for Z con¬ 
trol of the qubit are sent through a separate microwave 
line and combined with the offset current by means of a 
bias tee located at the base plate. The chosen attenua¬ 
tion in the Z pulse line sets the effective electron noise 
temperature in the center conductor to ~ 6mK above 
the base temperature. 


The quantum circuit is mounted in an aluminum box 
and is cooled down to about 20 mK in a dilution refrig¬ 
erator. Its sample box is enclosed in a cryoperm case 


II. PARAMETERS OF THE CONCENTRIC TRANSMON 
A. System Hamiltonian 


I Electronic mail: jochen.braumueller@kit.edu. 


The geometric inductance of the large loops in the con¬ 
centric transmon architecture forms a notable contribu- 






















































(a) 


according to 


2 



FIG. S2. Different schematic circuit diagrams of the con¬ 
centric transmon qubit. The black dot denotes the central 
qubit island, (a) Exact lumped element representation of the 
qubit circuit, (b) Simplified effective circuit diagram. The 
dc SQUID is treated as a single effective Josephson junction 
with critical current I^, connected in series with an effective 
inductance L and the qubit capacitance C. It is suggestive 
to identify the effective inductance with the inductance L of 
the full outer qubit ring, (c) Geometric representation of the 
circuit. The qubit capacitance is center symmetrically dis¬ 
persed. (d) Diagram motivating the choice of effective induc¬ 
tance L. Two half rings with inductance I//2, respectively, 
are connected in parallel, forming the geometric inductance 
L, = L/4. 


tion to the total qubit inductance. To account for this, 
we derive a generalized system Hamiltonian modeling our 
circuit. The lumped element representation of the qubit 
circuit is depicted in Fig. S2(a). The qubit capacitance is 
distributed center symmetrically into n small parallel ca¬ 
pacitances C/n and the total inductance Lg is split into 
a series of n inductors of inductance Lg/n. We approxi¬ 
mate this exact circuit by the effective simplified circuit 
diagram depicted in Fig. S2(b), consisting of a single ca¬ 
pacitor of capacitance C connected in series with an ef¬ 
fective inductance L and a single Josephson junction of 
effective critical current /*, valid in the regime L <C 
The total critical current of the SQUID is treated as an 
effective parameter I* = 1.5/c to account for model in¬ 
accuracies. The geometric inductance Lg we extract in 
this analysis is related to the effective inductance L via 
Lg = LjA. This can be motivated since Lg is formed by 
two half rings in parallel, see Fig. S2(c), (d). This choice 
of effective parameters is justified by a good agreement of 
the presented model with measured data and simulation 
results. In the following we present a detailed derivation 
of the system Hamiltonian based on the simplified circuit 
diagram depicted in Fig. S2(b). 

From current conservation at the two independent ac¬ 
tive nodes we get 

=/* sin02 = (SI) 

The phases (j)i in Fig. S2(b) correspond to the respective 
voltage drop across each of the components of the circuit 


fdt'U.it'). (S2) 

2e Jto 

The directed voltage drops in a closed network add up to 
zero according to Kirchhoff’s rule. In a closed loop, the 
integration constant in Eq. (S2) is the total magnetic flux 
penetrating the loop. Since the network in Fig. S2(d) is 
an open loop in the sense of the flux quantization law, 
we can write 


<('1 + 02 + <(>3 — 0- (S3) 

We eliminate 02 using Eq. (SI) and insert in Eq. (S3) 

03 -b arcsin(c03) = -0i, (S4) 

introducing the notation 

(S5) 

2e I*L ^ ‘ 

An approximate solution to Eq. (S4) can be written as 
1 . c 


03 = — 


(1 + c)' 


r0i - 


(l + c)= 


- sin 01 


(S 6 ) 


which is exact when 0i ^ 1 or c 1. Plugging into 
Eq. (SI) and writing 0i = 0, we reduce the set of equa¬ 
tions to a single equation of motion: 


Ar 0 = - — 

2e ^ 2eL(l + c)2' 


h 1 


2 eL (l-bc )2 
With the Euler-Lagrange equation 
d dC dC 


dt 


0 (^ 0 ) 
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we obtain the Lagrange function 
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Employing the definitions of the Josephson energy Ej = 
^/c, and the inductive energy El = (^) P^" 

rameter c becomes 


2e 1 

¥ 1.5A 
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El 

3Ej' 


(SIO) 


Using the conjugated variable for the charge number n = 
C ( 2 e )2 0^ arrive at the Hamiltonian 


6Ei 


H = AEc{n-ngY -Ej ^ 


(6Ej + 2Ei)2 
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9E2 


'(6Ej + 2El)2'^ 

with charging energy Eq = ^ j^C. 
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FIG. S3, (a) Fundamental qubit transition and lowest or¬ 

der multi-photon transitions. From the extracted transition 
frequencies woi/Stt = 7.6496 GHz, |LUo2/27r = 7.6094 GHz 
and Icuos/Stt = 7.5673 GHz we obtain three equations us¬ 
ing Eq. (S17). The transition peaks are measured at vary¬ 
ing excitation power (blue data points). Transparent data 
points correspond to a measurement of the full frequency 
range, (b) Frequency spectrum of the concentric transmon 
qubit with respect to the dc current applied to the flrrx bias 
line. The thick blue line denotes the measured qubit transi¬ 
tion frequency which is fitted (red) to the typical dc SQUID 
dispersion. The faint blue line is the two-photon transition 
to the second excited transmon level. Qubit measurements 
presented in the main text are taken at a;oi/27r = 6.85 GHz. 


B. Charging energy Ec, Josephson energy Ej and 
inductive energy El 


Expanding the constructed system Hamiltonian 


H = AEc{n —ngY — Ej cos ^ +, (S12) 


with 


E.j = El 


QEl 


{6Ej + 2ELy 


El = El 


9Ej 


i6Ej + 2ELy 


(S13) 


up to fourth order in (j) yields 

H = 4:Ec{n-ng)'^ + (^Ej/2 + Ei^ 0^- -|-const. 

(S14) 

We identify AEq = ^ and Ej/2 + Ei = to cast 

the Hamiltonian in Eq. (S14) into the standard form of 
the harmonic oscillator for terms of order 0{(jP). This 

TABLE SI. Extracted parameters of the concentric transmon 
qubit using data in Fig. S3. Errors in Ec and C are assumed, 
other errors are standard deviations as extracted from the 
fit. Ic denotes the total critical current of both Josephson 
junctions. The contribution Lg of the geometric inductance 
relative to the total qubit inductance Ltot at 6.835 GHz is 
Lg! Ltot = 10 %. 



Ej 

la 

Ec 

C 

El 


sweet spot 
6.85 GHz 

45 GHz 
29 GHz 

90 nA 
58 nA 

0.24 GHz 

81fF 

128 GHz 

0.64 nH 

A 

±12 GHz 

±23 nA 

±0.03 GHz 

±10fF 

±30 GHz 

±0.16 nH 


yields the phase operator 

^ = ^ («’''+a) (S15) 

\Ej/2 + Ei) 

in the representation of creation (annihilation) operators 
(a). Plugging into Eq. (S14) and taking into account 
the bosonic commutation relation, [a, d^] = 1, yields 

H = 4^Ec{Ej/2 + EL)a^a 

~ A{e7/2 + Ei) + (S16) 

and we can find the energies Eq^ of transmon levels \m) 
to be 

Eon = Ec{Ej/2 + EL)m- ~ (m^+m). 

(S17) 

Figure S3 (a) shows the fundamental qubit transition 
|0) -O' |1) and the lowest order multi-photon transi¬ 
tions 1/2 10) O 12) and 1/3 |0) e-i- |3), measured close 
to the current sweet spot of the spectrum, see arrow in 
Fig. S3(b). Equating Eq. (S17) for m = 1,2, 3 and us¬ 
ing measured transition data renders three equations for 
the transmon parameters Eq, Ej and Ei. In an iter¬ 
ative approach using this transition data and by fitting 
measured spectroscopy data, we can extract a suggestive 
set of parameters which are in good agreement with both 
data sets. The spectrum fit is carried out using Eq. (S17) 
evaluated for m = 1, see Fig. S3(b). Ej = Ej($) is a 
function of the applied magnetic flux $ with 
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^ 1 + (P tan^ ( TT -h A$ 
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(S18) 


The effective critical current /*° of both Josephson 
junctions is dependent on the external magnetic flux $, 
induced by the flux bias line. denotes the total crit¬ 
ical current of the concentric transmon at the flux sweet 
spot, $ = 0 and A$ is a phase offset. The square root 
term in Eq. (S18) accounts for a relative asymmetry be¬ 
tween the Josephson junctions®^, specified by the param¬ 
eter d = . A separate fit of the qubit spectrum 

indicates an asymmetry parameter of d = 0.32, causing a 
decrease in tunability range of the circuit. The minimum 
qubit frequency of the investigated sample is 6.3 GHz. 

Table SI summarizes the parameters extracted from 
the described fitting algorithm. The qubit capacitance 
deviates slightly from microwave simulation data of the 
geometry, predicting C = 58fF. Possible reasons for this 
discrepancy are an additional capacitive contribution by 
the coupling to the readout resonator in the experiment 
and the additional capacitive effect of the supply lines to 
the Josephson junctions which are not considered in the 
simulation. A geometric finite element simulation yields 
Lg ~ 0.2 nH, deviating by a factor of three from the fitted 
value. 

We want to point out that the presented model may be 
considered as a semi-quantitative approach to model the 
concentric transmon circuit. Within this framework, the 
assumptions made and the renormalization of parameters 
is justified by the good agreement of measured data and 
simulated data with the model. 


III. FINITE ELEMENT SIMULATION OF DEFECT LOSS 

To validate the microscopic analysis of defect loss con¬ 
tributing to qubit decay, we conduct a finite element sim¬ 
ulation of the concentric transmon geometry. We model 
the sample by assuming an oxide layer thickness of 3 nm, 


formed by AlOx (Sr = 11-5), surrounding the central 
island and the concentric ring electrode of the qubit. 
For the interface between substrate and aluminum we 
assume an effective dielectric constant of e^. = 6, ac¬ 
counting for the contribution of Si02 {Sr = 4) in the 
interface oxide. The simulation yields an electric field 
energy fraction of 2.8 • 10“^ residing in the oxide volume. 
With a phenomenological loss tangent®"^ '^AiOx = 3 • 10“^ 
of AlOx, this leads to a total effective loss tangent of 
Stlf = 8.4 • 10“^ attributed to dissipation by surface and 
interface oxide defects. Due to errors in JaiOx) the oxide 
layer thickness and Sr, this value is to be considered exact 
within a factor of two. 6tlf can be calculated to limit Ti 
at rj”TLF ~ ns which is in good agreement with the 
microscopically extracted value for a;g/27r = 6.85 GHz. 

In the same spirit we can extract a relative electric field 
energy ratio in the silicon substrate of 0.92. With an in¬ 
trinsic silicon loss tangent Ssi < 1 • 10“^, dissipation in 
the substrate can be neglected. This small value is sug¬ 
gested for instance by very high quality factors measured 
for TiN resonators on intrinsic silicon®^. 
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